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Introduction 

In  this  paper,  there  is  formulated  a  linear  theory  of  a 
three-dimensional,  elastic  continuum  which  has  some  of  the 
properties  of  a  crystal  lattice  as  a  result  of  the  inclusion, 
in  the  theory,  of  the  idea  of  the  unit  cell.  The  equations 
yield  wave-dispersion  relations  with  acoustic  and  optical 
branches  of  the  same  character  as  those  found  at  long  wave¬ 
lengths  in  crystal  lattice  theories  and  observed  in  neutron 
scattering  experiments.  Although  specific  solutions  are  not 
exhibited  in  detail,  it  is  apparent  from  the  form  of  the 
equations  that  there  will  be  interesting  surface  effects  under 
conditions  of  both  motion  and  equilibrium. 

The  unit  cell  may  also  be  interpreted  as  a  molecule  of 
a  polymer,  a  crystallite  of  a  polycrystal  or  a  grain  of  a 
granular  material.  The  mathematical  model  of  the  cell  is  a 
linear  version  of  Ericksen  and  Truesdell's  deformable  direc¬ 
tors  [1].  If  the  cell  is  made  rigid,  the  equations  reduce  to 
those  of  a  linear  Cosserat  continuum  [2]. 

The  method  of  derivation  of  the  equations  is  analogous 
to  one  used  in  deducing  two-dimensional  equations  of  high 
frequency  vibrations  of  plates  from  the  three-dimensional 
equations  of  classical  linear  elasticity.  By  the  same  tech- 


nique  as  that  employed  in  passing  from  high  order  theories  of 
plates  to  classical  theories  of  plates,  the  equations  are 
shown  to  reduce,  at  low  frequencies  and  very  long  wave-lengths 
in  isotropic  materials,  to  those  of  an  elastic  continuum  with 
potential  energy-density  dependent  on  strain  and  strain  gra¬ 
dient  and  kinetic  energy-density  dependent  on  velocity  and 
velocity  gradient. 

A  linear  form  of  Toupin's  generalization  [8,  Section  7] 
of  couple-stress  theory  [5-10]  is  obtained  by  eliminating  the 
difference  between  the  deformations  of  the  unit  cell  and  the 
surrounding  medium;  and  linear  couple-stress  theory  itself 
is  obtained  by  eliminating,  further,  the  symmetric  part  of 
the  strain  gradient.  Both  of  these  special  cases  are  also 
limited  to  low  frequencies  and  very  long  wave-lenths. 

1.  Kinematics 

Consider  a  material  volume  V,  bounded  by  a  surface 
S,  with  X^,  i  *  1,2,3,  rectangular  components  of  the 

material  position  vector,  measured  from  a  fixed  origin,  and 
x £  the  components,  in  the  same  rectangular  frame,  of  the 
spatial  position  vector.  The  components  of  displacement  of 
a  material  particle  are  defined  as 


ui  5  xi  - 


(1.1) 
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Embedded  in  each  material  particle  there  is  assumed  to  be  a 
micro-volume  V*  in  which  Xj  and  xj^  are  the  components 
of  the  material  and  spatial  position  vectors,  respectively, 
referred  to  axes  parallel  to  those  of  the  x^,  with  origin 
fixed  in  the  particle:  so  that  the  origin  of  the  coordinates 
xj^  moves  with  the  displacement  u.  A  micro-displacement 
jg'  is  defined:  with  components 

u!  2  x!  -  X!  .  '  (1.2) 

ill 


The  absolute  values  of  the  displacement-gradients  are 
assumed  to  be  small  in  comparison  with  unity: 


<<  1 


Chi? 

dX! 

1 


<  <  i  ; 


so  that  we  may  write 


(1.3) 


Uj  *  uj(xi»fc)'  (!•*) 

u«  «  uj(xi,x^,t),  (1.5) 

\riiere  t  is  the  time. 

Assume  that  the  micro-displacement  can  be  expressed 
as  a  sum  of  products  of  specified  functions  of  the  x^  and 
arbitrary  functions  of  the  x^  and  t.  As  an  approximation. 


du.j  du.i 

=  3JT  3  ai°j  > 


dul  du' 

axf  *  ^  5  *iuj  » 
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retain  only  a  single,  linear  term  of  the  series: 

uj  m  xk* kj  > 

where  is  a  function  of  the  and  t  only.  Then 

the  displacement-gradient  in  the  micro-medium  is 

Diu-  =  *ij  (i-7) 

i.e.t  the  micro-deformation  is  taken  to  be  homogeneous 

in  the  micro-medium  V1  and  non-homogeneous  in  the  macro- 
medium  V.  In  view  of  (1.3)2*  |^ij|  ^  7116  symmetric 

part  of  is  the  micro-strain? 

*(ij) 3  +  *ji>  t1-8’ 

and  the  antisymmetric  part  is  the  micro-rotation : 


*'Uj]  s¥*ij  -  *n) 


(1.9) 


An  alternative  interpretation  of  the  micro-deformation  is 
that  the  quantities  are  proportional  to  the  components 

of  the  displacements  of  the  tips  of  de formed) le  directors ,  as 
described  by  Ericksen  and  Truesdell  [l].  The  then 

are  the  components  of  rotation  of  the  Cosserat  triedre  [2,p. 122} . 

Ms  define  the  usual  strain  (now  the  macro-strain) 


§(Vj  +  djUi!'' 


(1.10) 


and  also  a  relative  deformation  (the  difference  between  the 
macro-displacement-gradient  and  the  micro-deformation) 


y.  .  =  d.u.  -  ii .  .  . 
Jij  13 


(1.11) 


and  a  micro-deformation  gradient  (the  macro-gradient  of  the 
micro-deformation) : 


Kijk  =  Vjk  • 


(1.12) 


All  three  of  the  tensors  e.  .,  y.  .  and  K.  ..  are  independent 

Xj  Xj  IJX 

of  the  micro-coordinates  x| .  Typical  components  of  7,^ 
and  are  illustrated  in  Pigs.  1  and  2. 

The  and  ^ ^  are  assumed  to  be  single  valued 

functions  of  the  x±,  leading  to  the  compatibility  equations 


emik®ntjaiajckt 

emijdiltjkt 

3i(ejk  +  ®jk  -  vjk) 


*ijk  ’ 


(1.13) 

(1.14) 


(1.15) 


where  to. .  is  the  macro-rotation: 

ij  - 


“ij  s  ?(aiuj  -  ajui> 


(1.16) 


and  is  the  alternating  tensor. 
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2.  Kinetic  and  potential  energies 


Let  the  micro-medium  be  a  parallelepiped  with  vol¬ 
ume  V1  and  edges  of  lengths  2d^  and  direction  cosines 
with  respect  to  the  axes  x! .  Let  x^  be  oblique  Cartesian 
coordinates  parallel  to  the  edges  d^,  respectively.  Then 

[3,  P.  153] 


x' 

i 


-  Vi 


9 


(2.1)x 


V'  a  8  !  I  -t.  .-t  |  |2  d  d  d  , 

•'  1]  lk"  1  2  3  3 


«•  -  1 1  .tlk  |  |*  dx"dx"dx“  . 

1.J  12  3 


(2.1), 

(2.D, 


Let  pM  be  the  mass  of  macro-material  per  unit 
macro-volume  and  let  p*  be  the  mass  of  micro-material  per 
unit  macro-volume.  We  define  a  kinetic  energy-density  (kinetic 
energy  per  unit  macro- volume) t 


P '  +  Uj )  («j  +  uj  )dv* , 


(2.2) 


where  the  dot  designates  differentiation  with  respect  to  time. 
Upon  substituting  (1.6)  and  (2.1)  in  (2.2)  and  performing  the 
integration,  we  find 

*  ■  4  <>Vj +  £ p  ,dkAj*tj  •  (s-i)l 


where 


P  5  Pm ■+  P'. 


(2.3). 


Kc  •'o=Vtki'u+  4D,4aA3,t3>=d?k  (2-3) 


p-  q 


peaks'*  Wu*k 


and  6ij  is  the  Kronecker  symbol.  In  the  case  of  a  cube  with 
edges  of  length  d  parallel  to  the  axes  of  x ! , 


'ij  =  4ij  •  d,  -  %  -  d,  3  d- 


(2.4) 


Then  the  second  term  in  (2.3) x  reduces  to  £  .  If 

the  material  is  composed  wholly  of  unit  cells,  =  0.  Then 


p»  -  p. 


For  the  potential  energy-density  (potential  energy 


per  unit  macro-volume) ,  we  assume  a  function,  W,  of  the 
forty-two  variables  cij »Kijics 


w  »  "Ui-j • 


(2.5) 


A  small,  rigid  rotation  of  the  deformed  body  is 
described  by  a  rotation,  *  constant,  of  the  macro-material 

and  an  equal  rotation  of  the  micro-material.  The  as¬ 

sociated  displacements  are 


uj  ■  *imij  •  uj  - 


(2.6) 
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The  addition  of  such  a  displacement  leaves  W  unchanged  since 


the  added  Y/_.  and  ^  are  zero. 


e  .  .  ,7  .  • 


The  assumptions  (1.6),  (2.2)  and  (2.5)  are  the  minimum 
that  will  lead  to  equations  which  yield  the  desired  dispersion 
relations  for  plane  waves:  including  longitudinal  and  transverse 
acoustic  and  optical  branches.  More  or  less  thant  (1.6),  (2.2) 
and  (2.5)  would  be  more  or  less  than  what  is  required. 

The  unit  cell  is  taken  to  be  a  parallelepiped  in 
order  to  represent  the  unit  cell  of  a  crystal  lattice.  How¬ 
ever,  another  shape  would  only  change  the  tensor  d^, .  Also, 
the  cell  can  be  interpreted  as  a  molecule  of  a  polymer,  a 
crystallite  of  a  polycrystal  or  a  grain  of  a  granular  material. 


3*  Variational  equation  of  motion 

We  write  Hamilton’s  principle  for  independent  varia¬ 
tions  6u^  and  between  fixed  limits  of  u^  and 

at  times  t  and  t  : 

o  i 

t  t 

6f  (^-7^)dt  +  /  b'W  dt  =  0  ,  (3.1) 

t„  1 

o  o 

where  &  and  "U T  are  the  total  kinetic  and  potential  energies: 

J  3  /  TdV,  TV=  /  WdV  (3.2) 

V  V 

and  bW  is  the  variation  of  the  work  done  by  external  forces. 

l 


wmmmmmmmesmmmmmmssm 


In  the  usual  way  [4,  p.  l66j,  we  find/  from  (2.3) 


and  (3.2)  , 


t  t  ! 

6f  7  dt  =  -/  dt  /(PU.5U.+  I  P'dj^tk^jk)dV.  (3-3) 
t/,  V 


As  for  the  variation  of  potential  energy,  we  first 


define 


_  _  aw 

°i3  •“  37- 


=  dW 

*^k=  '3*I* 


(3.4) 


(3-4). 


(3.4). 


Then 


W  -  Tij«€ij  +  0ij67ij  +  Uijk«Kijk, 


Tijai4Uj  +  °ij(ai5Uj  -  4*ij)  +  ^ijkW  * 
ai[(Tij  +  °ij)4ujl  -  °ij>4uj  -  aij4*ij 


+  ^ijk4*}*5  *  ^ijk^Jk 


Applying  the  divergence  theorem,  we  find 
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/  6WdV 
V 


•  "  -  i  di(Tij  +  °ij>tujdV  -  /('Vijk+-«jlc)4*’jkdv 

+  {  +  °ij>#ujdS  +  /  "^ijk^jkas  •  (3-5) 

The  form  of  (3*5)  is  the  motivation  for  the  adoption 
of  the  following  form  for  the  variation  of  work  done  by  ex¬ 
ternal  forces: 


45#'i“  l  Wv  + 1  tjk^jkdv  +  /  kj6ujds  +  f  Tjk5V 


dS  . 


(3-6) 


The  definitions  of  u.  and  V'-sv*  and  fche  fact  that 
the  integrands  of  the  volume  and  surface  integrals  represent 
variations  of  work  per  unit  volume  and  area,  yield  the  physical 
significances  of  the  coefficients  of  6u..  and  5^^*  Thus, 
f ^  is  the  body  force  per  unit  volume  and  t^  is  the  surface 
force  per  unit  area  (stress-vector  or  traction);  <1*^  is  to 
be  interpreted  as  a  double  force  per  unit  volume  [4,  p.  187 ] 
and  as  a  double  force  per  unit  area.  The  diagonal  terms 

of  and  are  double  forces  without  moment  and  the 

off-diagonal  terms  are  double  forces  with  moment.  The  antisym¬ 
metric  part  of  the  body  double  force  is  the  body 


couple.  The  antisymmetric  part  T 


[jl<] 


of  the  double  tractior 
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Tjk  is  Ith8  Cosserat  couple-stress  vector.  In  both  and 

the  first  subscript  gives  the  orientation  of  the  lever 
arm  between  the  forces  and  the  second  subscript  gives  the 
orientation  of  the  forces.  Across  a  surface  with  its  outward 
normal  in  the  positive  direction,  the  force  at  the  positive 
end  of  the  lever  arm  acts  in  the  positive  direction.  ("Positive” 
refers  to  the  positive  sense  of  the  coordinate  axis  parallel 
to  the  lever  arm  or  force).  Across  a  surface  with  its  out- 
word  normal  in  the  negative  direction,  the  directions  of  the 
forces  are  reversed. 

Substituting  (3.3),  (3-5)  and  (3.6)  in  (3.1),  and 
dropping  the  integration  with  respect  to  time,  we  obtain  the 
variational  equation  of  motion: 


f(dixij  +  ^i°ij  +  f j  _  pUj )6UjdV 


+  £<Vijk+  °jk+  V  3  p'dj/tk>8Viv 


+  /  ttrn.(T.j+aij)]«ujdS  +  /(Tjk-ninijk)«*jkdS-° 


(3.7) 

4.  Stress-equations  of  motion  and  boundary  conditions 

From  the  variational  equation  of  motion,  there  fol¬ 


low  immediately  the  twelve  stress-equations  of  motion: 
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di<*ij  +  °ij>  +  £j  "  P“j  - 


(4.1) 


Vijk  +  °jk  +  *jk  -  3  p'dt/tk-  ^-l); 


and  the  twelve  traction  boundary  conditions: 


t.  =  ni(T..  +  0ij), 


Tjk=  "i^ijk  ' 


(4.2) 


(4.2) 


In  view  of  (3-4),  (4.2)  and  the  significance  of 

e..,  y..,  k. ...  appropriate  terminology  appears  to  be  Cauchy 
13  1]  XJK 

stress  for  T^j»  relative  stress  for  and  double  stress  for 

The  twenty-seven  components  of  are  interpreted 

as  double  forces  per  unit  area.  The  first  subscript  of  a 
designates  the  normal  to  the  surface  across  which  the  component 
acts;  the  second  and  third  subscripts  have  the  same  significance 
as  the  two  subscripts  of  T^*  Typical  components  of 
are  illustrated  in  Fic.  2. 

The  linear  equations  of  a  Cosserat  continuum  [2]  are 
obtained  by  setting  ^(j_j)  *  0*  Th®11  °(ij)  *  Tij  ^i(jk)*^ 

and  there  remain  (the  Cosserat  couple-stress)  and  a j ± j  j  s 

which  has  been  regarded  as  the  antisymmetric  past  of  an  asym¬ 
metric  stress  However,  in  the  present  theory,  the  Cauchy 

stress,  Tij>  is  symmetric  and  °[ij]  i®  4he  antisymmetric 
part  of  an  asymmetric  relative  stress  j . 
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Beside3  containing  the  linear  equations  of  a  Cosserat 
continuum  as  a  special  case,  Eqs.  (2.5),  (3.4)  and  (4.1)  also 
include,  as  low  frequency,  very  long  wave  length  approximations, 
linear  versions  of  the  equations  of  couple-stress  theory  [1-10] 
and  Toupin's  generalization  of  couple-stress  theory  (8,  Section  j), 
These  are  considered  in  Sections  9-12. 

If  additional  terms  were  retained  in  the  series  ex¬ 


pansion  (1.6)  of  the  micro-displacement  uj,  higher  order 


stresses  would  appear.  In  addition  to  stresses  corresponding 
to  double  forces  per  unit  area,  there  would  be  stresses  cor¬ 
responding  to  n-tuple  forces  per  unit  area.  All  of  the  latter 


would  be  self  equilibrating;  whereas,  of  the  only  the 


Li(j*) 


are  self  equilibrating. 


5.  Constitutive  equations 


For  the  potential  energy-density  we  take  a  homoge¬ 


neous,  quadratic  function  of  the  forty-two  variables 


* 


ijk' 


W  "  2°ijkt€ijckt  +  2  bijkt7ij7kt  +  2  aijktmnKijkKtmn 


+  d 


ijktra7ijKktm  +  fijkfcn*ijkefcn  +  gijkt7ij€kt*  t5’1) 


Only  I  X  42  X  43  “  903  of  the  42  X  42  ■  1764  coefficients 
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in  (5.1)  are  independent.  The  number  of  coefficients,  the 
relations  among  them  and  the  number  of  independent  one3  are 
given  in  the  following  table: 


cijkt  “  ck'ij  "  cjikl  : 

9  x  9  -  60 

=  21 

bijk*  =  hktij  : 

9  >9  -  36 

=  45 

a  .  — »  a  • 

ljktmn  Imnijk  * 

27  x  27  -  351 

=  378 

dijk?m  : 

9  x  27 

=  243 

£  a  £  • 

ijklm  ijkmt 

9  x  27  -  81 

=  162 

gi jk l  =  gij<k  s 

9  x  9  -  27 

-  54 

903 

From  (5.1)  and  ( 3 . 4 ) : 


T  s*  c  .  .t  .  .  + 

.  pq  pqao  13 

pq  ypqi3  13 

11  =  f  .  .c, .+ 

^-pqr  pqri}  ij 


g .  .  7  .  .  + 

yi3pqri3 


b.  .  7 .  .  + 

i3pq'i3 


dijpqrYij+ 


fijkpqKijk  * 
dpqijkKijk  , 
apqrijkKijk* 


(5.3) a 

(5.3) , 

(5.3) 3 


In  the  case  o£  a  centrosymmetric,  isotropic  material 
(referred  to  as  isotropic  in  the  sequel)  the  number  of  in¬ 
dependent  coefficients  is  greatly  reduced.  As  there  are  no 
isotropic  tensors  of  odd  rank,  an<*  ^ijktm  must  vania 
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The  remaining  coefficients  must  be  homogeneous,  linear  func¬ 
tions  of  products  of  Kronecker  deltas.  There  are  three  inde¬ 
pendent  products  of  two  Kronecker  deltas  and  fifteen  independent 
products  of  three  Kronecker  deltas.  Hence 


cijkt  =  *6ij6kt  +  Ui6ik6jt  +  tt*.fiit6jk  '  (5-4) x 


b .  ,  s  b  4.  .i,  .+  b  4..  A.,  +  b  4.  ,4.,  ,  (5.4) 

13k!  i  i]  kl  2  lk  ]t  a  xi  ]k  ’  w  '2 

q.  ..  t  *  g  6  .  .6.  ,+  g  6  .,6  . ,  +  g  6  .  .6  ..  ,  (5-4) 

yi3kt  13  kl  ■* 2  ik  ]t  ■* 3  it  3k  *  3 


a..,  .  8  a  4.  .4.  ,6  +  a  4.  .4,  6  +  a  6  .  .  6,  6 , 

ijktmn  i  13  kt  mn  2  13  km  nl  3  13  kn  tm 


+a  6  6 .  .6  +  a  6., 6.  6  ,+  a  6  6 .  6, 

4  3k  it  mn  s  3k  ini  nt  s  3k  m  tm 


+a7akiSjiamn+  aBakiajmant+  aBaki6jnain> 


+a  6  . .  6  .  6,  +  a  6. .6.  6.  +  a  4.  .6.4. 
10  it  jin  kn  11  3t  km  in  12  kt  im  3n 


+a  6., 6.  6.+  a  6..6V  6.  +  a  6V.6.  6._ 
^3  iv  jn  Km  14  3 ^  Kn  xm  is  K%*  in  3® 


(5.4) 


The  conditions  (5.2)  require  the  six  relations 


u.  ■  p.  =  it  ,  g  *  g  > 

*1  2  2  3 


s  ■  s  i  s  ■  a  i  3  ■  n  >  b  *  3  t 

1  0  2  9  3  T  11  12 

leaving  eighteen  independent  coefficients.  Thus,  the  potential 
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onorgy  density  reduces  to 

W  =  —  .  .c  •  .  +  uc.  ^  b  7.  .7,  ,  +  ■!  b  7,  .7,  . 

2  11  33  i3  X3  2  x'li-  3  j  2  2rij  rij 

+ 1  Vij’ji +  5 jj  +  9,  <'ij  +  Ve« 

*  ^iikVjj-"  a2ltiik>‘jkj+  I  ViikKjjk+  I  VijjKikk 

+  a3lcijjKkik+  2  aaKijiKkjk+  2  a1„ltijkKijk+  a11KijkKjki 

+  2  ai3KijkKikj+  2  an'tijk,tjik+  4  ai=Kijk'tkji'  (5-5) 
and  the  constitutive  equations  become 


x  3  X5  e .  .  +  2ac  +  g  6  y  ■  •  +  g  (y  +  y  ) , 
pq  pq  ll  pq  yx  pq'll  y«>wpq  |rqp' » 


(5.6) 


°pq  "  ^x6pqCii+  2g=Sq+  blWii  +  Vpq  +  Vqp  *  <5*6>S 

^pqr  **  ai^Kiip6gr+  xrii6pq^+  a2^xiig6pr+  Kiri6pq^+  a3xiir6pq 

+  a  k_..6 _ +  a  (k  .  .6  +ie.  .6  )+  a  k.  .  6  +  a  k 

4*^x1  qr  5V  qxx  pr  xpi  gr'  a  lqx  pr  i0  pqr 

+  a  +  x___)+  a  k  +  a,  k__  +  a  x„„^.(5.6) 

ix  rpq  qrp'  x3  prq  X4  qpr  is  rqp  : 


6.  Displacement-equations  of  motion 


We  may  obtain  twelve  equations  of  motion  on  the 
twelve  variables  and  by  first  inserting  (1.10), 


>:v*>  ^jrjusr><\f4 £'- 
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i 


(l.ll)  and  (1.12)  into  the  constitutive  equations  and  then 
the  latter  into  the  stress-equations  of  motion  (4.1).  There 
is  no  necessity  to  assume  spatially  homogeneous  material  pro¬ 
perties.  In  fact,  the  coefficients  of  elasticity  and  the  den¬ 
sities  mav  be  taken  as  periodic  functions,  of  oblique  coordi¬ 
nates  parallel  to  the  edges  of  the  unit  cell,  of  periods  2d^. 
This  would  represent  the  periodic  structure  of  a  crystal 
lattice.  However,  the  equations  would  then  be  highly  intrac¬ 
table;  whereas  some  of  their  main  features  arc  exhibited  with 
macro-homogeneous  material  properties,  at  least  for  wave¬ 
lengths  greater  than  the  dimensions  of  the  unit  cell.  The 
isotropic  case  is  especially  simple;  but  nevertheless  it  still 
contains  many  of  the  novel  properties  of  the  macro-homogeneous 
material. 

In  the  case  of  isotropy,  the  constitutive  Eqs.  (5.6) 
apply  and  also  (2.4).  Then,  for  the  macro-homogeneous,  iso¬ 
tropic  material  the  equations  on  u^  and  are 

(n  +  2gg+  b2)  +  (x  +  lJL  +  2gi+  b^+  b3)^juj 

-(gA+  bA)^jj-(g2+  b2)djV'ji-(9s+  b3)dj^ij+  fi~  P"i*  (6*1* 


+  aioWij+  9A“k6ij+  92(di“ 

+  "Av  <'kk>6ii+  Wr  ^ij)+  b3(ajui- 


.+  b  .u. ) 
3  3 


*3*'  13 


rd’*ii 

(6.2) 


7-  Micro-vibrations 


Consider  solutions  of  (6.1)  and  (6.2)  of  the  form 


ui  -  fi  “  *«  - 0  •  *a  -  Aije 


rat 


(7.1) 


where  the  A^j  are  constants.  Then  (6.1)  are  satisfied 
identically  and  (6.2)  become 


bi6ijAkk  +  bzAij  +  b3Aji  “  3  P'd2a)2Aij  »'  (7*2) 

which  admit  the  following  solutions: 
dilatational  mode: 

A»  -  ’  \.  ’  Aij  "  °  *  1  *  1  ’  (7'3^ 

<o*  -  3(3^  +  ba  +  b^/p'd*  ,  (7.3)„ 
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(7.4) 1 

(7.4) 2 

(7.5) 2 
(7-5)2 

(7.6) 1 

(7.6) , 


The  restriction  of  the  potential  energy-density  to 
be  positive  definite  requires 

3b  +  b  +  b  >  0  ,  b  +  b  >0,  b  -  b  >  0  . 

X  2  3  23  '23 

Bence  cn^,  a>g  and  o>r  are  real  frequencies.  The  corresponding 
modes  are  analogous  to  the  simple  thickness-modes  of  vibration 
of  a  plate.  Just  as  the  latter  are  independent  of  the  coordi¬ 
nates  in  the  plane  of  the  plate,  so  are  the  micro-modes  inde- 
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pendent  of  the  coordinates  x^  of  the  three-dimensional  con¬ 
tinuum  with.  micro-3tructure.  Extensional  and  flexural  waves 
in  a  plate  couple  with  thickness-modes,  at  high  frequencies, 
to  form  the  highe..  ^ranches  of  the  dispersion  relations  for 
a  plate.  Analogously,  we  may  expect  longitudinal  and  trans¬ 
verse  acoustic  waves,  in  the  three-dimensional  continuum  with 
micro-structure,  to  couple  with  the  micro-modes  to  form  optical 
branches. 

8.  Plane  waves,  long  wave-length 

Consider  solutions  of  (6.1)  and  (6.2)  with  and 
zero  and 

ui  =  ui(xx»fc)»  *ij  =  ^ij(x1»t)  •  (8*1) 

By  means  of  linear  combinations,  the  twelve  equations  may  be 
composed  into  three  independent  equations  and  three  independent 
systems  of  three  equations  each: 


“3P'dSW(8-2) 

as  (8.2)  except  that  i® 


;; 

w 


shear  optical  I; 


$ 


?3 

H 

’V 

U 


(a  +  a  )d  d  if,  \  - 

'  lO  Z3/  1  1  (23/ 


shear  optical  II:  the  same 
replaced  with  if  -if  : 

*  22  33 
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rotational  optical; 


(a  —  a  )b  d  i!/  —  (b  —  b  )i/  —  ■hp,d2^/  t  (8.3) 

V  lO  13;  lY[23]  1  2  3)y[23]  3 H  ^[23]'  ' 


longitudinal  system : 


k  d  d  u  -  k  d  $  -k  d^  =  pii  , 
li  1  1  1  12  1  11  1'  1  1 


k  d  u  +k  3  3/  -  k>  +  k  33^=  -^p'd2^0  ,  )  (8.4) 

SI  1  1  22  1  i’ll  22  11  23  1  1  2r  T  11*  ( 


k  du  +k  3  3^  +  k  3  3  ^  -  k1  i  =  p  *  d2tf 

31  1  1  32  1  lrll  33  1  1T  33 


1  D 

where  «  -r  and 

3  rii  *  rn  rn 


k  «  X  +  2p.  +  2g  +  4g  +  b  +b  +b  , 

11  ^  *1  *2  1  2  3 

k  »2a+a+a+-§a  +2a  +  ■§  a  +  a  +  a  ,  . 

22  2  3  8  2  lO  11  2  13  14  15 

k  »  6a  +  2a  +  a  +  9a  +  6a  +  a  +  3a  +  2a,  +  3a  +  a,  +  a,  » 

33  1  2  3  4  5  8  lO  11  13  14  15 

k  ■  k  =*3a+2a+a+3a+a+  2a, ,  +  a,  +  a  , 

23  32  1  2  3  5  8  11  14  15 


k  ■  k  «  3g  +  2g  +  3b  +  b  +  b  , 

Sl  13  1  32  J  1  2  s' 


2g  +  b  +  b  , 

2  2  3 


\  +  bJ' 


3  (31>l  +  ba  +  ba), 


transverse  syatem  It 


‘u‘A*. '  k»A*(l2f  ki.Vtuf  PV 

\AV  k«W(xa)-  k^(x2)+  S2.SiV[  >2]‘  3p,d2?(x2))( 


k  au  +  k  a  a  \+ k  a  a  *  r  *'■ ^p*d  *n.2i 

31  1  2  32  1  1  (12)  33  1  1  [12J  33  [ 12 J  3  l 12 i 


where 


11  +  2g  +  b  , 
^  2  2 


2a  +  a  +  a  +  2a^  +  2a  +  2a  +  a  +  a  , 

2  3  8  lO  11  13  14  15 

-  2a  +  a  +  a  +  2a.  -  2a,  -  2a  +  a  +  a  , 

2  3  8  lO  11  13  14  15 


k  *  a  -  a  +  a  -  a  , 

32  3  B  14  15 


k  ■  k  ■  b  -  b  , 

31  13  2  3 


k  *  k  -  2g  +  b  +  b  , 
12  21  2  2  3 

k'  «  2(b  +  b  )  , 

22  .  V  2  3 


k'  *  2(b  -  b  )  , 

33  2-3 


transverse  system  XX ;  the  same  as  transverse  system  I  except 
that  V  *.l2)  and  *(l2j  are  replaced  with  *(l3) 

and  if,  ,,  respectively. 
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If,  now,  (8.1)  are  specialized  to  the  plane  waves 

UjL  -  Aiiexp[i(4xi-  cut)],  expUfex^-  cut)],  (8.6) 

four  dispersion  relations  (<u  vs.  §)  result: 
shear  optical  waves  ( SO ) ( twice ) : 


•i  p •  d2ui2  »  b  +  b  +  (a  +  a  )£2  ?  (8*7) 


2  3  *  XO  13 


rotational  optical  waves  (RO) 
—  n  •  d2m2  a  h  — 


f  p'd*m2  -  'b  -  ba  *  (aio  -  ai3U 


longitudinal  waves  (LA,  10,  LDO) : 


(8.8) 


$2  -  pen2 

k  t 

12’ 

k  £ 

13^ 

k  £ 

21 

k  52  +  k'  -  i  p»d2u>2 
22  22  « 

k  *2 

23 

k  * 

3X 

k  42 

k  $2  +  k* 

32 

33  33 

(8.9) 


transverse  waves  (TA;  TO,  TRO)  (twice): 

1 1 22  n/i\2  ip  S 


k  -  par 
xi* 


k  4 
ax 


k  i 
ist9 


^aX? 


+  K*-  3  p’dS 


k  ( 

13 


is  e* 

88 


k„«*  +  |p’« 


(8.10) 
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These  dispersion  relations  are  similar  to  those  en¬ 
countered  in  a  second  order  theory  of  extensional  waves  in 
plates  [11].  The  relations  (8.7)  and  (8.8)  are  like  that  for 
the  second  face-shear  mode  in  the  plate.  The  relations  (8.9) 
and  (8.10)  are  like  that  for  the  coupled  extensional,  thickness- 
stretch  and  symmetric  thickness-shear  modes  in  the  plate;  i.e., 
one  acoustic  and  two  optical  branches  in  each  of  (8.9)  and  (8.10). 

In  the  dispersion  relations  (8.7)  and  (8.8)  for  the 
non-coupled  modes,  there  are  cut-off  frequencies  and 
respectively,  at  which  the  group  velocity  (du/d£)  is  zero. 
Positive  definiteness  of  W  requires 

\o*\.>0*  aiO-aX3>0-  t8’11* 


Hence,  the  frequencies  increase,  from  cut-off,  with  increasing 
real  wave  numbers.  Below  the  cut-off  frequencies,  the  wave 
numbers  are  pure  imaginary  with  cut-off  values 


(8.12) 


respectively,  at  zero  frequency. 

The  behavior  of  the  acoustic  branches  in  (8.9)  and  (8.10),  at 
low  frequencies,  is  described  by  oj^,  ca”,  to”1  (i  *  1  for  longi¬ 
tudinal  and  i  -  2  for  transverse):  the  values,  at  oi  »  0  and 
£  a  0,  of  the  first,  second  and  third  derivatives  of  cd  with 
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respect  to  <; .  We  find 

=  Vf  ,  =  °  ,  c>»»  =  3v.(t^  -  h2),  (8.13) 

where 

v*  =  (X  +  2£L)/p  ,  v2  =a/p  ,  (8.14) 

89p  (39+  2g2)2  # 

X  +  2p.  =  X  +  2vi-  2(b^+  bg)  "  3(3bi+  b2+  ba),(8,15^i 

29* 

V-  "  ^  “  b" . +  b .  >  (8.15)2 

23 


2(a  +  a  +  a  +  a  +  a  )/(X  +  2£),  I2 

X  2  3  4  u  if 


2(a3+a4)/JI,(8.l6) 


h2  »  p»d2[ai2  +(ot  +  6)2)/3P,  h2  -  p»d2(l+82K6p, (8.17) 


a 


bj(39i+  29aA 
3bi+  v  bay 


2g2 

p  -  1+ftV> 


(8.18) 


ax-  |[(l+0)(3a48)ai+(l+2ap+02)a2-  ^(l40)(l-2o-P)a3-(l-P)(3DH0)a5 

-  ~(l-8)(l+2a40)aB+  2apaii-a(l-p)ax4+a(l40)ai5]  , 
a2-  |{-(l-aa-p)(3a40)ai-  |[l-(2a40  )2]aa+  |(l-2j-3)2a3+(3cH£)2a4 
+(3a-t0)(l+2a40)a8+  |(l+aa40)2aa+  a(3a+2p)axo+2x(a+£)axx 
+  a(3a+2fJ)axa+  a(Ha-tf  Ja^-  a(l-a-p)axs). 
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3a=  il-(l-P3)a2+  i(l-^)aa3+  , 

a4=|[(l-«2)aio-(l-e2)(aii+ai3)+|(14P)zai<+4(l-p)2ai5], 

a3=  -|[-(l-32)aio+(l+332)aii+(l-«2)ai3-  |( l+2@-302)ai4-  J. 

Positive  definiteness  of  W  requires  £L,  "X  +  2p[,  X?  >  0,  while 
h?  >  0  by  inspection . 

It  may  be  seen  that  the  limiting  group  velocities 
v^  are  less  than  those  that  would  be  calculated  from  the 
strain-stiffnesses  X  +  2\i  and  p..  Thi3  phenomenon  is  due  to  the 
compliance  of  the  unit  cell  and  has  been  found  in  a  theory  of 
crystal  lattices  by  Gazis  and  Wallis  [12].  Inasmuch  as  (u"  =*  0, 
the  group  velocities  at  zero  frequency  are  maxima  or  minima. 

Which  one  occurs  depends  on  whether  o>|M  is  greater  or  less 
than  zero;  and  this,  in  turn,  depends  on  whether  X2  is  greater 
or  less  than  h2.  Now  and  h2  are  positive  quantities 

that  are  length-properties  of  the  material  —  depending  on 
stiffness  ratios,  density  ratio  and  the  size  of  the  unit  cell. 
Although  d  is  probably  smaller  than  the  X^,  the  density 
ratio  p'/p  and  the  stiffness  ratios  a  and  0  can  make 
either  X  or  h  the  greater.  Hence,  as  the  frequency  in¬ 
creases  from  zero,  both  group  velocities  can  increase  or  both 
can  decrease  or  one  can  increase  and  the  other  decrease  —  depend¬ 
ing  on  the  properties  of  the  material.  There  is  no  analogue 
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in  the  theory  of  homogeneous  platoa  because  they  do  not  have 
multiple  stiffnesses  and  densities.  The  phenomenon  does  occur, 
however,  in  sandwich  plates  [13]  and  it  has  also  been  found 
in  a  theory  of  crystal  lattices,  with  complex  interatomic  in¬ 
teractions,  by  Gazis  and  Wallis  [14]. 

At  the  short  wave-length  limit  (£  **  “),  the  asymp¬ 
totic  values  of  the  group  velocities  of  the  acoustic  branches 
are 

+  k  -[(k  -k  „ ) 2 

22  33  *  22  33' 

from  (8.9)  and  the  same  expression,  with  k  replaced  by  k 

9 

from  (8.10).  These  can  be  much  smaller  than  (8 . 13) ^  if 

4(k  k  -  k2  )  «  (k  +  k  )2  ,  (8.21) 

22  33  23'  22  33 7 

4(k  k  -  k2  )  «  (k  +  k  )2  ,  (8.22) 

22  33  23  22  33' 

which  appear  to  be  possible. 

Regarding  the  optical  branches  in  the  longitudinal 
and  transverse  systems,  the  former  have  long  wave-length  cut¬ 
off  frequencies  u>8  and  <ud  while  the  corresponding  quan¬ 
tities  for  the  latter  are  u>8  and  o>r.  Thus,  the  two  systems 

have  one  cut-off  frequency  in  common.  One  of  the  two  modes 


4k=3r 


1 


(8.20) 
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is  shear  and  the  other  is  equivoluminal  extension.  As  in  the 
case  of  plates,  the  group  velocities  of  all  four  optical  bran¬ 
ches  are  zero  at  the  long  wave-length  cut-off  frequencies  ex¬ 
cept  in  the  unlikely  circumstance  of  coincidence  of  cut-off 
frequencies  within  a  system  (b^  =*  0  for  longitudinal  modes 
or  b  53  0  for  transverse  modes).  Another  parallel  to  the 

3 

situation  in  plates  is  that,  as  *  increases  from  zero,  the 
behavior  of  the  optical  branches  is  very  sensitive  to  small 
changes  in  the  ratios  of  material  properties.  One  possibility 
is  that  both  lower  optical  branches  have  phase  and  group  veloc¬ 
ities  of  opposite  sign;  i.e.,  diminishing  a,  with  increasing 
£ .  With  further  increase  of  £ ,  the  absolute  values  of  the 
group  velocities  would  pass  through  maxima,  then  drop  to  zero 
and  then  increase;  i.e.,  the  dispersion  curves  first  would 
have  a  point  of  inflexion  and  then  a  minimum. 

A  sketch  of  a  possible  configuration  of  the  real 
segments  of  the  dispersion  curves  is  shown  in  Fig.  3*  T^e  four 
lowest  branches  (TA,  LA,  LO,  TO)  are  remarkably  similar  to 
those  obtained  by  Brockhouse  and  Iyengar  [15,  Fig.  5]  from 
measurements  of  neutron  scattering  in  germanium. 

At  the  corresponding  stage  in  the  development  of 
equations  of  high  frequency  vibrations  of  plates,  it  is  ex¬ 
pedient  to  introduce  correction  factors  to  compensate,  as  well 
as  possible  within  the  framework  of  the  theory,  for  errors 
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introduced  by  the  restrictive  assumption  regarding  the  varia¬ 
tion  of  displacement  through  the  thickness  of  the  plate.  The 
analogous  restriction,  here,  is  the  assumed  homogeneous  de¬ 
formation  of  the  unit  cell.  The  values  of  the  correction  fac¬ 
tors  are  obtained,  in  the  theory  of  plates,  by  matching  appro¬ 
priate  points,  slopes  and  curvatures  of  the  dispersion  curves 
with  the  corresponding  quantities  obtained  from  an  exact  so¬ 
lution  of  the  three-dimensional  equations.  Since  the  analogue 
of  the  latter  does  not  exist,  in  the  present  case,  such  an 
adjustment  cannot  be  made.  An  alternative  is  to  use  experimen¬ 
tal  data. 

There  is  another  aspect  of  the  theory  of  plates  that 
should  be  mentioned.  It  is  possible  for  a  thickness-mode 
with  n  +  m  nodal  planes  to  have  a  frequency  lower  than  one 
with  n  nodal  planes.  For  example,  in  an  isotropic  plate, 
the  thickness-shear  mode  with  two  nodal  planes  has  a  frequency 
lower  than  that  of  the  thickness-stretch  mode  with  one  nodal 
plane  if  Poisson’s  ratio  is  greater  than  one-third.  Thus,  a. 
better  approximation  is  obtained  if  a  sufficient  number  of 
terms  is  retained,  in.  the  series  expansion  of  the  displacement, 
to  accomodate  this  contingency  [11],  The  analogue,  in  the 
present  case,  is  the  possibility  of  the  appearance  of  a  micro¬ 
mode  with  frequency  lower  than  that  of  the  dilatational  micro¬ 
mode  (7.3)  if  additional  terms  are  retained  in  the  series 
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expansion  (l.6).  However,  because  of  the  complications  that 
would  ensue,  such  a  step  does  not  appear  to  be  warranted  at 
this  time. 

9.  Low  frequency,  very  long  wave-length  approximation! 

Form  I 

This  section  and  the  following  three  are  devoted  to 
discussions  and  derivations  of  equations  of  motion  simpler 
than  (6.1)  and  (6.2)  but  limited,  in  application,  to  much 
lower  frequencies  and  much  longer  wave-lengths. 

As  noted  previously,  when  thickness- shear  and  thick¬ 
ness-stretch  deformations  and  the  associated  inertias  are 
taken  into  account  in  the  theory  of  plates  [16, 17,18),  thickness 
nodes  of  vibration,  analogous  to  the  micro-modes,  are  obtained 
as  well  as  flexural  and  extensional  modes  analogous  to  the 
transverse  and  longitudinal  acoustic  modes.  At  low  frequencies, 
in  comparison  with  the  frequencies  of  the  thickness-modes,  and 
at  long  wave-lengths,  in  comparison  with  the  thickness  of  the 
plate,  the  coupling  of  the  flexural  and  extensional  modes  with 
the  thickness-modes  is  negligible.  As  the  frequencies  of  the 
flexural  and  extensional  modes  approach  zero,  the  thickness- 
shear  deformation  approaches  aero  but  the  thickness-stretch 
deformation  does  not;  rather ,  it  is  the  stress  associated  with 
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thickness-strotch  that  approaches  zero.  Thus,  the  antisymmetric 
and  symmetric  parts  of  deformation  and  stress  have  to  be  treated  dif¬ 
ferently  in  passing  from  high  frequency  equations  to  the  class¬ 
ical,  low  frequency  equations.  To  obtain  equations  valid  at 
low  frequencies  in  the  case  of  flexure,  the  thickness-shear 
deformation  is  made  to  approach  zero  by  passing  to  a  limit  as 
the  associated  modulus  of  elasticity  approaches  infinity. 

The  product  of  the  two  is  indeterminate  and  this  leaves  the 
thickness- shear  stress  indeterminate  in  the  constitutive  equa¬ 
tions.  In  the  case  of  extension,  the  thickness-stress  is  set 
equal  to  zero  and  the  resulting  constitutive  equation  is  used 
to  eliminate  the  thickness-strain  from  the  remaining  equations. 

In  both  flexure  and  extension  of  homogeneous  plates, 
the  thickness  velocities  are  set  equal  to  zero  in  the  kinetic 
energy,  for  the  low  frequency  approximation,  because  their 
contributions  are  negligibly  small  at  the  low  frequencies  to 
which  the  resulting  equations  are  restricted  owing  to  the 
suppression  of  the  thickness-shear  deformation  and  the  omission 
of  the  thickness-stretch  stress  [18,  19] •  The  same  is  not 
true  of  non-homogeneous  plates.  For  example,  in  a  sandwich 
plate  the  rotatory  inertia  of  the  facings,  about  the  middle 
plane  of  the  plate,  can  be  of  paramount  importance,  even  at 
low  frequencies,  for  certain  combinations  of  stiffness  ratios, 
density  ratio  and  distance  between  facings  [13 I. 
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Now,  the  thickness  velocities  are  analogous  to  the 
micro-velocities  j »  thickness  of  the.  plate  i3  ana¬ 
logous  to  dimensions  of  the  unit  cell;  the  thickness- 

shear  deformation  is  analogous  to  the  antisymmetric  part  of 
the  relative  deformation  (^[ij]^*  *die  thickness-shear  moduli 
are  analogous  to  the  and  fc^e  stress  associated  with 

thickness-stretch  is  analogous  to  the  symmetric  part  of  the 
relative  stress  (a(ij))*  T*ie  known  process  of  descending  from 
high  frequency  equations  of  plates  to  low  frequency,  long 
wave-length  approximations  can  serve  as  a  guide  to  the  treat¬ 
ment  of  analogous  terms  in  the  equations  of  the  elastic  con¬ 
tinuum  with  micro-structure.  However,  regardless  of  the  process 


or  of  the  theory  of  plates,  the  test  of  the  validity  of  the 
resulting  equations  of  motion  is  that  they  yield  the  same  dis¬ 
persion  relations  that  are  found  in  the  limit  as  o  ■*  0,  £  "*  0 
for  the  acoustic  branches  of  the  general  equations.  Thus,  in 
the  isotropic  case,  the  values  of  ru?,  tuV  and  ta!n ,  in  (8.13)» 
must  be  reproduced  exactly.  Attention  will  be  confined,  here, 
to  this  case  because  it  is  much  simpler  than  the  anisotropic 
form  and  because  (8,13)  are  available  for  the  final  test. 

Inspection  of  (8.13)  shows  that,  as  in  the  case  of 
sandwich  plates,  it  is  not  permissible  to  discard  the  micro- 
velocities  Their  effect  is  contained  in  the  ha  and, 

as  remarked  in  the  discussion  following  (8.19)*  ha  can  be 
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leas  than  or  greater  than  tz  depending  on  stiffness  ratios, 
the  density  ratio  and  the  dimensions  of  the  unit  cell.  Hence, 
omission  of  the  $ ^  would  preclude  the  reproduction  of  the 
low  frequency  behavior.  The  remainder  of  the  process,  however, 
can  follow  the  analogy  with  homogeneous  plates. 

Thus,  we  let 

°(ij)  =  °’  O'1) 

ba  '  ba  "*  ”  ’  ’  0  '  (9-2) 

and  proceed  to  find  the  effect  of  these  assumptions  on  the 
remaining  terms. 

The  isotropic,  constitutive  equations  for  T  and 
a^,  separated  into  symmetric  and  antisymmetric  parts,  are 


t  »  X6  € . .  +  2u€  +g5  y . .  +  2g  y  r  \  ,  (9*3) 

pq  pq  u  ^  pq  p<r  ii  2  (pq)  *  v 

°(pq)’  9i4pq€ii+  2Vpq+  Vp-^ii  +  <V  b3h(pq)  '  {9A) 

°[pq]*  (V  b3hIpq]  •  (9’5) 

Then,  with  (9.2),  °[pq]  is  indeterminate  in  (9.5)  and,  with 

(9.1) ,  (9.4)  may  be  solved  for  7(pq)  in  terms  of 

■»(M)  M  +  (l_e)£M’ 

where  a  and  P  are  given  in  (8.18). 


(9.6) 
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With  regard  to  c>  we  note  first  that,  since 


and  7 


(pq3 


7  =  3  U  -  tit 

*pq  p  q  *pq 


is  now  zero,  we  are  left  with 


^[pq]  “pq  '  ^(pq)  "  epq  "  7(pq);  (9-7) 


or,  using  the  expression  (9.6)  for  7(pq)  >  we  have 


'(pq) 


=  ai  e .  .  +  . 

pq  n  pq 


(9.8) 


Accordingly,  *ijk  =  =  d.*(jk)  +  ^[jk]  reduces  to 

Kijk  aKUt6jk  +  ■g(l4^)Kijk“  2(1_^Kikj'  ^9‘9^ 


where 


*ijk  ~  3iSj“k  =  Kjik* 


(9.10) 


Thus,  the  part  of  the  potential  energy-density  that  is  a  func¬ 
tion  of  *^jk  becomes  a  function  of  rijks  the  second  gradient 
of  the  displacement.  The  eighteen  components  of  “i^jk  may 
be  resolved,  in  more  than  one  way,  into  tensors  whose  components 
are  independent  linear  combinations  of  the  80  bbat 

other  forms  of  the  energy-donsity,  for  the  low  frequency  ap¬ 
proximation,  are  possible.  These  are  treated  in  Sections  11 
and  12 . 
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Upon  inserting  (9*9)>  (9-6)  and  (9*2)  in  (5*5)>  we 
find  Form  I  of  the  low  frequency  approximation  for  the  poten¬ 
tial  energy-density: 


W->  W  =  2  \eiie^  +  +  a2KijjKikk 


+  a3KiikKjjk  +  \KijkKijk  +  asKijkKkji  '  (9-H) 


where  A  and  £  are  given  in  (8.15)  and  a^...as  in  (8.19). 
The  appearance  of  these  coefficients  is  preliminary  evidence 
of  the  validity  of  the  process. 

We  define  new  stresses: 
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The  variational  equation  of  motion  is  now  obtained 
from  Hamilton's  principle  with  independent  variations  6u^ 
alone  since,  by  (9*7)1  and  (9-8)»  the  ip ^  are  no  longer 
independent  of  the  u^. 

The  variation  of  the  potential  energy-density  is 


6W  -  + 


-  W  SijAYv 

-  ^jt^jk-  ai*ijk>4uk1~VV4i5iJk>4V5iGiJkV“k)- 

(9.14) 


Hence 


/awav./nj  (T  jk-diSi  jk)  Kds-/a  j  ^ik-^i^ijk )  6ukdv+£ni>Ii  jkd  j6ukds 

(9.15) 

NOw,  in  the  last  integral  of  (9»15)>  the  variation  is 

not  independent  of  6u ^  on  S:  only  its  normal  component 
njdjfiu^  is  independent.  We  separate  the  latter: 

ni?ijkSj4uk=  ni5ljkDj*V  n#ijknjD5“k«  <9-l6> 


where 


Di  3  (sjt-  njnt)aV  D  2  "A  • 


(9.17) 


The  terms  in  (9.16)  may  be  resolved,  further,  in  more  than  one 
way.  Xn  this  section  we  follow  Toupin  [8]  and  reserve  an 


alternative  resolution  [9]  for  Section  12.  Thus,  for  the  first 
term  on  the  right  hand  side  of  (9.16),  which  contains  the  non- 


independent  variation 


we  write 


'VijkV'V  Dj<niiIijk6uk)-niDj'lijk5uk-(Djni)iIijk6uk-  (9-lB> 


The  last  two  terms  in  (9.18)  now  contain  the  independent  varia¬ 
tion  6u^.  For  the  preceding  term,  we  note  that,  on  the  surface 
S, 


Dj(n^ijk6uk>=(Dtnl>njni*1ijk6uk+nqeqpmap(emtjniniiIijk6ukM9.19) 


By  Stokes's  theorem,  the  integral,  over  a  smooth  surface,  of 

the  last  term  in  (9*19)  vanishes.  If  the  surface  has  an  edge 

C,  formed  by  the  intersection  of  two  portions,  S  and  S  , 

1  2 

of  S,  Stokes's  theorem  gives 


/nqeqpmyemljntnl?ijkSuk)dS  *  *  Jnimj‘LijkJ6ukds>  (9'20> 

s  c 

where  m^  *  e^j  s^n^  and  the  sm  are  the  components  of  the 

unit  vector  tangent  to  C.  The  bold  face  brackets  [  ]  in 

(9;20)  indicate  that  the  enclosed  quantity  is  the  difference 

between  the  values  on  S  and  S  . 

1  2 

Finally,  we  note  that,  in  the  first  surface  integral 
in-(9.15)»  we  may  write 

njai?ijk  "  njDJ?ljk+  ninjD^ijk' 


(9.21) 


Then,  assembling  the  results  in  (9-15)-(9.21)  we  find 


/  WdV 
V 


'  iaj(?jk-  di5ijk)6ukdv 
+  /[nj?jk-ninjDSijk-2njD.ii..k+(ninjDtnt-Djni 

+  {ninjiIijkD4ukdS  +  i  lnimjPijkl4ukdS- 


Wijk>4ukds 

(9.22) 


This  form  suggests,  for  the  variation  of  work  done  by  external 
forces. 


*  fFk6ukdV  +  -^k61^133  +  /\D6ukds  +  ^E^Su^ds.  (9.23) 

V  s  s  c 


As  for  the  kinetic  energy,  the  micro-velocity 
in  (2.3) must  be  replaced  by  a  linear  function  of  macro-veloc¬ 
ity  gradients: 

-  hijkAit  •  o-2") 

where 

hijkt  H  5^4ikSj(.'4ii4jk)4a4ij4k(.+  ?“(aik4jf.+4it4;jk^  ’ 

(9.25) 

so  as  to  satisfy  (9.7)  and  (9*8).  Then  the  kinetic  energy- 
density  (2.3)  becomes 
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■?pVj+s(IWMv 

■  2  P.VJ  +i  3p[p'apkinn(a">'1n)'5lcl-  §  V  P'^fa^VV 

(9.26) 


where 


^2  2  j2  v  i*  _  j2 

pkmn  “  jt n£qpk  jqmn  “  mnpk 


I  a!(Vta-Vk.t2a(5!,2s>W!(VkntV‘k.)l 


(9.27) 


The  total  kinetic  energy  is 


T-  /Tdv=/ [ipu^-  ^p(p'd=lOTnam.-.n)0k]dv+/  IP'npd^Cd^l^dS, 


ln.„  3= 


V  V 


from  which 


(9.28) 


6/>dt-  dt/ipi^-  i  8p(p'a;tam»man)]s'Vdv 

t_  t  V 


-/dt  /  fp1npa*tann(Dman+nmDun)eukds  . 


(9.29) 


The  variational  equation  of  motion  is  formed  from  (9.22),  (9*23) 
and  (9.29),  from  which  follow  the  stress-equations  of  motion 
and  boundary  conditions! 
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V  PV  3  V^pkamW  (9-3°)  j 


Vjk-ni"jBiIijk-2njDiiiijk+(nin3Dtn<.-Djni^ijk 


+  3  P'Vptan^mV-  "mD“n>=  V  O'30)* 

ninjPijk=  \<  (9-3°)3 

Inin^ijkJ“V  ^•30). 


The  displacement-equations  of  motion  are  obtained 
by  first  replacing  2e^j  with  djUj  +  ^jui  and  *^ijk  with 
Si^jU^,  in  (9.13),  and  then  substituting  the  latter  in  (9*3°) x> 
The  result  is 

( A+2u )  ( l-t^V2 )  W- u-£  ( 1-1*  V2 )  VXVXU  +  F  *  p(u-h^W.u+  h^XVXu)  , 


(9.31) 

where  the  t ?  and  h|  are  defined  in  (8.16)  and  (8.17). 

Omitting  the  body  force  and  taking  the  divergence 
and  curl  of  (9-31),  we  find  the  equations  governing  the  prop- 
agation  of  dilatation  and  rotation: 


-(l-h^jv.a  ,  (9  32Ja 

v^i-X^)  7=  yxu  .(i-W)v*5  ,  (9.32), 

w*lere  tha  vf  ara  dafinad  In  (8.14).  For  tha  plana  wavs s 
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(^•Jtf>Vxu)=(A,A)exp[i(^n.r  -  cot)],  (9-33) 

the  dispersion  relations  are 

-  v?52(l  +  t?52)/( 1  +  h?e2),  (9.34) 

from  which  follow  exactly  the  properties  (8.13).  Thus  the 
validity  of  the  approximate  equations  for  low  frequencies  and 
very  long  wave  lengths  is  established.  The  dispersion  rela¬ 
tions  (9.3*0  are  illustrated  in  Pig.  3  by  the  dashed  curves 
labelled  L  and  T. 

10.  Relation  to  Toupin* 8  generalization  of  couple-stress 
theory 


The  theory  of  elasticity  with  couple-stresses,  which 
is  considered  in  [5-10],  is  based  on  the  same  kinematics  as  is 
classical  elasticity;  but  the  potential  energy-density  is  as¬ 
sumed  to  be  a  function  of  the  strain  and  the  curl  of  the  strain 
instead  of  the  strain  alone.  In  the  linear  theory,  the  compo¬ 
nents  of  the  curl  of  the  strain  are  the  same  as  the  components 
of  the  gradient  of  the  rotation:  eight  independent  linear 
combinations  of  the  eighteen  components  of  the  second  gradient 
of  the  displacement.  For  the  equilibrium  case,  Toupin 

[8,  Section  7l  has  generalized  the  theory  to  include  all  eighteen 
components.  If  the  inertia  terms  are  omitted,  (9*30)  ere  iden- 
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tical,  in  form,  with  Toupin's  Eqa.  (7.8)-(7. 11)*.  However, 

*  . 

although  the  form  is . the  same,  these  are  some  significant 
differences.  Equations  (9.3°)  pertain  to  a  low  frequency, 
very  long  wave-length  approximation  to  the  equations  of  a 
material  with  micro-structure  and  the  effect  of  the  micro¬ 
structure  survives,  in  both  the  potential  and  kinetic  energy- 
densities,  through  the  contribution  of  the  symmetric  part  of 
the  relative  deformation.  This  part,  can,  in  fact  be 

traced  to  T.  .  and  p...v,  in  (9*30).  through  the  coefficients 
a  and  p  in  W.  Similarly  the  contribution  of  to 

the  acceleration  terms  in  (9*30)  can  be  traced  through  the 
coefficients  a  and  P  in  h^^  and  then  in  0n 

the  other  hand,  Toupin's  equations  do  not  stem  from  considera¬ 
tions  of  microstructure. 

The  equations  of  the  material  with  micro-structure 
can  be  reduced  to  those  of  a  material  without  microstructure 


(i.e.,  to  a  micro-homogeneous  material)  by  causing  the  micro¬ 
medium  to  merge  with  the  macro-medium.  This  may  be  accomplished 
(just  as  readily,  in  this  case,  for  the  anisotropic  as  for  the 
isotropic  medium)  by  passing  to  the  limit  as 

^ijk“t  °°  *  "^ij  +  ®  ’  (iO.l) 

*  Note  that,  by  definition,  Toupin's  is  symmetric  in 

the  second  and  third  indices,  whereas  ftijv  is  symmetric  in 
the  first  two  indices.  Note,  also,  thatJ  b  should  be  replaced 
by  -b  in  Toupin's  Eqs.  A,B,  (7*9)  (7.19)* 
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in  stead  of  (9.1)  and  (9.2);  and,  at  the  same  time,  p  ->  p*, 

* 

so  as  to  remove  the  distinction  between  micro-and  macro-den¬ 
sities  (i.e.,  pM  ■*>  0  in  (2.3)^).  Then 


*  ai“j  ’  diai“k  3  3ijk  > 


(10.2) 


instead  of  (9*7) (9-8)  and  (9*9).  Accordingly,  from  (5.1), 

_  1®  X  Q  M  **  o  M 

W  W  -  2  cijktGijckt+  2  aijktmnKijkKtmn+  fijktmKijkGtm(l03) 


and  we  define  new  stresses 


=  dw  ~o 
Tij  -  ■  Tji » 


(10.4) 


jo  s  j«l.  „  JO 
^ijk"  a*ijk 

Also,  in  the  kinetic  energy-density  (9.26), 


(10.4). 


^pkmn  dpm^kn' 


(10.5) 


where  the  d^ft  are  again  given  by  (2.3) 

In  the  isotropic  case  d^L  **  d26pjn  and  W°,  T°j>  V'-ijk 

have  the  same  form  as  W,  T^j,  but  the  coefficients  are 

A,  a,  a  ...a  instead  of  T,u,  a  ...a  . 
r  1  S  X  s 

The  formulation  of  the  variational  equation  of  motion 
proceeds  as  before  and  we  arrive  at  the  stress-equations  of 
motion  and  boundary  conditions} 
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di>Iijk)+  V  P'V  3  y>,aym"k)'  (10-6), 

"j’jk-  niniD'^ijk-  2njDiP°jk+(ninjVt-  Djni)?ijk 

t3p,”pdyD.VnAl,',k  ’  (10-6)2 

ninjP°jk  *  *£  -  (10-6)3 

InimjPijkJ  *  =k  ’  (10-6>4 

Without  the  acceleration  terms,  (10.6)  are  now  pre¬ 
cisely  the  linear  form  of  Toupin's  equations;  and  (9. 30), 
without  the  acceleration  terms,  differ  only  in  that  the  co¬ 
efficients  in  ISr  and  W°  (and  hence  in  the  stresses)  have 
different  meanings.  However,  with  the  acceleration  terms  in¬ 
cluded,  the  forma  on  actually  different  i  the  fourth  rank  tensor 
cl^kran  in  (9-30)  is  replaced  by  the  second  rank  tensor  d^  in 
(10.6)  so  that  there  are  fewer  coefficients  in  the  latter. 

In  the  isotropic  case,  it  is  only  necessary  to  let 


X  X,  u  -*  p.,  o  ■*  0,  0  -*•  1,  (or  g^-*  0, g^-* 0)  ( 10 . 7) 

in  the  equations  of  Section  9»  to  reach  the  equations  of  the 
micro-homogeneous  medium.  Thus,  (9»3i)  reduces  to 

(X+2u)(l-?*7a)W.u-^(l-?|7a)y>^>qj  +  P»  p*(l-  ^d8^)*?,  (10.8) 


-46- 


where  l 2  and  -t2  are  obtained  from  t2  and  i2,  in  (8.l6), 
by  employing  (10.7).  It  will  be  observed  that  the  left  hand 
sides  of  (9*3l)  and  (10.8)  have  the  same  form,  but  the  right 
hand  sides  have  different  forms  because,  with  (10.7),  h2=h2=  ^d2. 

12  3 

The  dispersion  relations  for  plane  waves,  from  (10.8)  are 


>f  -  v||=(i+?|5a)/(i+i  an=). 


(10.9) 


where  v2  =*  (X+2p.)/p‘,  v2  =*  p./p*;  and  the  properties  at  u>  =  0, 
£  =*  0  are 


»  vA,  u>£  =  0,  <u"'  =  3Vi(?|-  -|d2)  (10.10) 


instead  of  (8.I3). 

The  difference  between  the  equations  of  the  micro- 
homogeneous  medium  and  the  equations  of  the  low  frequency, 
very  long  wave-length  approximation  is  similar  to  the  difference 
between  plane  strain  and  plane  stress;  or,  more  appropriately  in 
the  present  context,  to  the  difference  between  equations  of  low 
frequency  extensional  vibrations  of  plates  with  the  thickness 
of. the  plate  constrained  and  not  constrained  to  remain  constant. 
In  the  case  of  equilibrium,  the  difference  is  solely  in  the 
physical  interpretation  of  the  elastic  stiffnesses.  With  stiff¬ 
nesses  determined  by  experiments  falling  within  the  restrictions 
of  the  equations,  the  two  equilibrium  theories  would  be  indis¬ 
tinguishable.  For  example,  the  numerical  quantity  that  would 
be  assigned  to  the  stiffness  £L,  in  one  theory,  would  be 
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assigned  to  y.  in  the  other.  In  the  case  of  motion,  however, 
there  is  an  essential  difference  as  the  number  of  coefficients 
is  not  the  same  in  the  two  theories  unless  g^  **  g^  =0,  which 
is  analogous  to  zero  Poisson’s  ratio. 


11 .  Low  frequency,  very  long  wave-length  approximations 
Form  II 


As  noted  before,  the  eighteen  components  ^chu^ 
may  be  arranged  in  independent  linear  combinations  which  form 
tensors.  One  such,  indicated  by  Toupin  [8,  p.  404],  is  the 
gradient  of  the  strain: 

4ijk  3  ■  ^aiaj“k  -  aAuj>  -  4ikj  •  (ua)- 

The  potential  energy-density,  for  the  low  frequency, 
very  long  wave-length  approximation,  may  be  expressed  as  a 
function  of  and  by  setting 

*ijk  -  4ijk +  Vi  -  Ka  <u-2) 

in  (9.11),  with  the  result: 


W->'W  -  I  ue..e..+  34 


2  Atii6jj+  aeijeij+  ai*iik*kjj+  aa*ijj*ilck+  a3Kiik*jjk 


A  A  A  A  A  A  .  _  \ 

+  a4KijJcKijk+  a5Kijk*kji  »  (13**3) 


where 


a  3  2a  -4a  ,  a  =»  -  a  +  a  +  a. 

1  1  3  2  X  2  3* 

a  =*  4a  ,  a  =  3a  -  a  ,3  *  -23  +  23  .  (11.4) 

3  3  4  4  5  5  4  5  \  / 


New  stresses  are  defined  by 


-  _  dW  j 

ij  be ^ j  =  ji  ' 


(11.5), 


^ijk“  “  ^ikj 

°Kijk 


whence 


(11.5). 


?pq  "  ^apq€ii  +  ^pq  ' 


“  -k  a  (6  k  .  .+  20..  +  4k..)  +  2a  4  JL, 

pqr  2  x'  pq  rxx  qr  xxp  rp  qxx'  2  qr  gpj 


(11.6) 


+  a  (6  k..  +  6  k..  )  +  ?.a  r  +a  (k  +k  J.(ll.6) 

3V  pq  Hr  priiq'  4 ^pqr  sv  rpq  grpr  v  *  ' 


The  variation  of  the  potential  energy-density  is 


6W  -  ?ijic6«ijk*  *ij^i6uj+  fijk^idj6uk 

"V^jk*'ai^ijk^6ukl"Sj^jk"ai^ijk^6uk  +  ai^ijkSjfiuk^  (H-7) 


Now,  (11.7)  has  the  same  form  as  (9.1*0  and  the  kinetic  energy 
a  a 

density  (9.26)  is  unchanged.  Hence,  the  variational  equation 


.  sbimt  .V  <B3.ovw-m.«w««»  im- -t 


-49- 

of  motion  has  the  same  form  as  in  Section  9  and  leads  to  bound¬ 
ary  conditions  like  those  in  (9*30)  and  stress-equations  of 
motion: 

W-  dAjk>+  Fx=  "v  3  •  (u-8) 

Recalling  that  =  £jik  wherea3  =£ikj>  the  quantity 

in  parentheses  on  the  left  hand  side  of  (9*30)  x  is  n°t  sym¬ 
metric  but  the  corresponding  quantity  in  (11.8)  is  symmetric. 

The  latter  is  a  more  convenient  form  for  the  introduction  of 
a  stress  function  of  the  Airy  type. 

To  get  the  displacement-equation  of  motion,  sub¬ 
stitute  (11. 1)  and  (1.10)  in  (11.6)  and  the  latter  in  (11.8). 
The  result  is 

(x+giXi-t^v2 )w-u-£( l-t^v2 ) vxTxu  +  p  =  p(u-h^w.u  +  h|y^»)  , 

(11.9) 

where 

t*=2(3  +3  +3  +3  +3  )/(X+2j),  t2=(3  +2a  +3  )/£fl.  (11.10) 

X  2  3  4  5  2  3  4  5 

In  view  of  (11.4),  2?  *  *l|,  so  that  the  displacement-equations 

of  motion  (11*9)  and  (9.31)  are  identical. 

12.  Low  frequency,  very  long  wave-length  approximation: 
Form  III 

For  some  purposes  it  is  advantageus  to  separate  the 
curl  of  the  strain  (or.  the  gradient  of  the  rotation) : 

& 

r 
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r,  .  =  e.,  b,e  .  =  4  e.,  d.d.u  . 
ij  j£m  t  mi  2  j-tm  l  t  m' 


(12.1) 


from  33  K^j  is  the  part  of  that  gives  rise 

to  couple-stresses.  The  double  stress  is  separated,  thereby, 
in  to  a  non-self-equilibrating  part  and  a  sel-f-equilibrating 
part.  Now,  =  0  .  Hence  Ic^  has  only  eight  independ¬ 

ent  components.  (They  are  the  components  of  the  dyadic  k 
in  Reference  [9 1 - )  The  remaining  ten  linear  combinations  of 
the  d.d.u.  were  considered  separately  by  Jaramillo  [20]. 

i  3  K 

They  can  be  expressed  as 

Kijks  Ki.jk+  3  °i'.j*kt.+  3eitkKjt'  3't)i‘ijuk+  S]'\ui) • 

(12,2) 

Thus>  Sijk  “  “kij  ’  *jki  ”  !jik!  ie-  *ijk  is  fully  oy"- 
metric. 

The  potential  energy-density,  for  the  low  frequency, 
very  long  wave-length  approximation,  may  be  expressed  as  a 
function  of  ic^  and  by  setting 

Kijk  "  *ijk  ’  3  eUj*kt”  3eitkKjt  (12-3) 


in  (11.3).  The  result  is 

W  ^  W  "  2  *eii€j j+  ?€ijeij+  2diKijKij+  2d2KijKji 

+  2  axKiij“kkj+  *2Kijk*ijk+  ?eijk*ij*kU  »  (12*4) 
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where 

l8d  =  -  2a  +  4a  +  a  +  6a  -  3a 

1  X  2  3  4  5 

iea2  -  23i  -  432  -  S3,  3si  =  2(3i  +  r  +  Zj, 

5  3  a  +  a  ,  3f  =  a  +  4a  -  2a  . 

2  4  5  X  2  3 

(12.5) 

The  definitions 

-  _  dW  - 

T  .  .  =  *r -  3  T  .  .  . 

13  de^ 

(12.6)x 

-• 

_  aw  A 

>Hj  =  "  °- 

(12.6)_ 

2 

t' 

^ijk=  “fkij’M.jki  “^jik> 

0Kijk 

(I2.6)a 

vr 

where 

is  the  couple-stress  deviator,  lead  to 

a. 

t  *  A6 _ +  2Ue^  , 

pq  pq  pq  pq 

( 12.7) A 

“l 

% 

*  43  k__  +  43  k__  +  fe„.K. . 

^pq  x  pq  2  qp  pqx  xji 

• 

(12-7)a 

V] 

% 

ft 

h 

^pqr" 

“i<*iir4pq+  *iip4qr+  Jiiq4rp)  +  252*pqr 

V* 

1 

+  3  Kij(Spqeijr+  4qreijp+  4rpeijq>  • 

(12.7)a 

i 

The  variation  of  the  potential  energy-density  now 

r1 

takes  the 

form 
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6W  =  Sytiy*  Sijk**ijk. 


djt^jk-3i5ijk!6uk)-dj<Ti3-3iiIIjk)4Vai(iIijkauk!'(12-8) 


where 


“*  =  i 


^ijk=  2  ejk-^it+  ^ijk 


(12.9) 


Again,  (12.8)  has  the  same  form  as  (9-14)  and  so  we  can  find 
stress-equations  of  motion  and  boundary  conditions  of  the 
same  form  as  (9-30)  *  but  with  replaced  by  *ijk‘  Such 

a  form  of  the  boundary  conditions  cannot  be  compared  directly 
with  the  results  of  Reference  (9]  because,  there,  one  of  the 
independent  variations  was  taken  to  be  the  tangential  com¬ 
ponent  of  rotation  -  which  here  is  embedded  in  the  normal 
derivative  D6u^.  To  get  the  alternative  form  of  the  boundary 
conditions,  we  can  return  to  (9*16)  and  further  resolve: 

D6\  «  26wjnjeijk+  Dk(ni6ui)-(Dkn.)6u.+  n^e^  (12.10) 

where  w^(=  eitm^^um^  is  the  rotation  and  enn  (not  summed) 
is  the  normal  component  of  strain  ninj€j_j  •  Then  integration 
by  parts  and  application  of  the  divergence  theorem  and  Stokes's 
theorem  leads,  in  the  notation  of  Reference  [9]*  to 
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6w  »  -/(y. r  +  \  yxv-d-  v-u-vj-iudv 

y  ^  M  *•  M  ' 


+/^n-T+  |nx  (V*u-^\Inn)-( V*jl)  •n-n*7x[nx(n*^+n ^ *nn)  ]}  -6udS 
S 


+/[n*ilxn+2nx(n*ji.n)xn]*(6wxn)ds  +  /nn  :j£ .  nfie  ds 
s  ~  ~  ~  -  ~  ~  ~  ~  ~  s~~  ~  1111 

+H|  i^nnS  +(sx  n)  -(n •£  +  n-u -nn)J -Suds  , 


(12.11) 


where  s  is  a  unit  vector  tangent  to  the  edge  C. 

The  variation  of  work  done  by  external  forces  is 
now  taken  to  be 

6  j?  -  /  F-5udV +/(P*6u+Q-6wx  n+Rde _ )dS  +  £  E-6uds  .  (12.12) 

i  ^  ~  -  ~  ~  nn  C  *"  ~ 

Here  Q  is  the  tangential  component  of  the  couple-stress  vector 
and  R  is  a  double  force  per  unit  area,  without  moment,  normal 


to  S, 


The  variational  equation  of  motion  then  yields  the 


stress-equation  of  motion 


v-T  +  yxy.^-y-jj.y  +  f  «pu-  -iv.(p'd2:Vu)  (12.13) 


and  the  boundary  conditions 


fl *5 +-§ a x(V •£-$!  n ) - ( V .5 )  •£-£} •VxfpxCp -jj+o  •  J •  nn )  ]+  -p'n-J2 :VU-P, 


(12.14) A 
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(12.14) 2 

(12.14) 

3 

(12. 14) ^ 

The  displacement  equation  of  motion  is  obtained  by 
substituting  (1.10),  (12.1)^  and  (12. 2)^  in  (12.7)  and  the 
latter  in  (12.13).  The  result  is 

(*+2u)(l-£2V2)VV-u-u(l^2V2  )  VXVxu+F=p  (u-h2W .  u+h^VxVXu ) ,  (12.15) 

Where 

*M3a  +  2a  )/(X+2H),  t®=( 3S  +a +2a  -f)/^L.  (12. 16) 

*  *  2  2  X  X  2 

In  view  of  (12.5)  and  (11.4), 

Z?  *  t|  =*  l?  »  ,  say;  (12.17) 

so  that  all  three  forms  of  the  low  frequency,  very  long  wave¬ 
length  approximation  yield  the  same  displacement-equation  of 
motion.  Necessary  and  sufficient  conditions  for  positive 
definiteness  of  W(»W»W)  are 

II  >  0,  3X  +  2JL  >  0  , 

a  >  o,  -a  <  a  <  a  , 

i  x  a  N  x  * 

a  >  0,  3a  +  2a  >  0,  5  <  0  . 

8  18 

Bence  4.®  >  0  and  we  have  already  seen  that  h|  >  0  . 


n  .{Ixn+2nx(n  -jjl  •  n)x  n  =  Q, 


nn:£-n  =  R. 

t+m  *0 


it  ^nnS  +{s^Mn*5+ n-J-nn)J  ,  3 


2  . 


(12.18) 


The  equations  that  were  considered  in  Reference  [9] 
are  obtained  from  those  of  this  section  by  setting  k,  =»  0 

1  jJC 

in  the  potential  energy-density  and  p»  =  0  in  the  kinetic 
energy  density.  Accordingly,  the  limitation  to  low  frequencies, 
long  wave-lengths  and  large  dimensions  is  more  severe  than 
was  apparent  previously. 


13.  Solution  of  the  approximate  equations  of  equilibrium 


In  this  Section  it  is  proved  that  any  solution,  u, 
of  the  equation 


( A+2u ) ( 1-t2 V2 ) W • u-u ( 1-t2  V2 ) VXVXU  +  F  =  0  , 

(13.1) 

in  a  region  V  bounded  by  a  surface  S,  can  be  expressed  as 

u 

-  B  -  |(kr^V2)V[r.(l-t=V2)^  +  Bq], 

(13.2), 

where 

H(l-t2V2)V2B  =  -  F  , 

(13-2)a 

ftCl-t2?2)^  =  r-Cl-t^jF  -  U2V-F, 

(13-2)3 

=  (X+u)/(x+2£L) 

and  r 

is  the  position  vector. 

Consider  a  field  point  P(x,y,z)  and  a  source  point 
And  define 
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whs  rs 

“(x-O2  +  (y-n)2  +(z-C)2  ,  dVQ=  d£dT}d£ . 
Then  7 aU  -  u  [3,  p.210],  or 


77-U  -  VxVxU  =»  u 

MM  M  M  M  M  M 

(13.^) 

Define 

r(t  =  7-U  ,  H  *  -  7  XU,  (7*H  -  0). 

(13.5) 

Then,  from  ( 13  -  4 ) , 

u  =  7 *  +  7  X  H  ,  7-H  ■  0 

~  r  mm'  ~  M 

(13.6) 

which  is  Helmholtz's  resolution.  Substituting  (13.6)  in 
(13.2)  ,  we  have 

jlT^kU-iJV8)^  +(l-t272)7xH]  +  P  -  0,  (13.7) 

where  k  ■  (X  +  2£L)/JL  .  Define 

•“3T  /  -t 

2  /  r"xe  x  2[k(l-t2V2)^+{l-t272)yXH]QdVQ  (13.8). 
Then  [3,  p.  210] 

(l-*,272)B»  -k(l-t272)7^  +(l-t=Va)VxH  (13.9) 
and,  from  (13.9)  and  (13.7)> 

vL(l--t2Vz)V2B  -  -  F  . 

Also,  the  divergence  of  (13*9)  yields 


(13.10) 


(13-11) 


( 1-t2 V2 )  V  •  B » =  k(  l-l2 V2 )  72^ 

Define  2k^*  =  r .(l-t272)B»  .  (13-12) 

Then,  using  (13.10) ,  we  have 

2kZ(l-l*V2)V21r*=  U^V*F-  r.(l-t2V2)F+2vL(l-t2V2)7-B‘ .  (13-13) 


Define 

BQ  =  2k[i/-i>*) 

(13.1^) 

and  find 

H(l-t^Va)V2Bo  =  r-(l-l2V2)F  -  U27-F 

(13.15) 

by  using 

(i3.ll)  and  (13.13).  Also,  from  (13- 1*0 

and  (13.12), 

2k^»r  * ( l-t272 ) B ' +  Bq  . 

(13.16) 

Now,  define 

p  =  p'-'t^7y’B'-k(l-t2V2 )Vtj, 

(13.17) 

By  (13.II),  V.D  -  0;  which  is  a  necessary  and  sufficient 
condition  for  the  existence  of  a  function  H*  such  that 
VX  H*  -  D,  i.e., 

V X  H*=»  B'  — t2  W*B ' -k(  1— 1272 ) 'Tjf  .  (13*18) 

mm  mm  '  m  2^"*  **  1  ~ 

Making  use  of  (13*9)  an(*  (13«H)  we  fin^ 

(l-t2V2)7XH*«(l-t*V2)7>«  .  (13.19) 

B?’h  7X  H-VxH* 


Next,  define 


(13*20) 
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and  note  that,  by  (I3.I9)  and  (13.20), 

(1-^V2)B"  *  0,  7-B"  *  0.  (13.21) 

From  (13.20),  using  (13.I8)  and  then  (13.16), 

VxH=.B,"+B'-i=W.B'-  |(l-t=Vs)7[r.(l-t=72B')+  Bj.  (l3.22) 


Then  substitute  (13.22)  and  (13.16)  in  (13.6^  to  get 

A(ki-t*7*)7(r.(l-t|7*)B-+B0J.  (13.23) 


Finally,  define 


B  =  B'  +  B" 

M  M  * 


(13.24) 


In  view  of  (13.21),  we  may  write  (I3.23)  and  (13.10)  in  the 
form  of  (13. 2)^  and  (13.2)^?  and  (13.15)  is  already  in  the 
form  (13.2)  .  Thus  Eqs.  (I3.2)  are  a  complete  solution  of 
(13*1).  If  ^  “0*  (13.1)  reduces  to  the  equilibrium  equation 
of  couple-stress  theory  and  (13.2)  is  the  solution  found  in 
[9].  If  both  l *  and  l*  are  zero,  (13.1)  is  the  classical 
equation  of  equilibrium  with  body  force  and  (13.2)  is  the 
solution  found  in  [21].  If,  in  addition,  the  body  force  is 
zero,  (13.2)  is  Papkovitch's  solution  [22],  The  proof  follows, 
generally,  that  in  [23]  but  with  an  improvement  as  a  result 
of  an  illuminating  criticism  by  E.  Sternberg. 
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1 4 .  -Concentrated  force  according  to  the  approxima t e 
equations 

In  an  infinite  region  V,  let  the  body  force  be 
zero  outride  a  finite  region  VQ  which  contains  the  origin 
and  a  non-vanishing  field  of  parallel  forces  P.  A  concentrated 
forrT  is  defined  by 

P  =  lim  /  F  dV  ,  (U.l) 

v  -*  0  V  ■  w 
o 

In  [9]  it  was  shown  that,  in  an  infinite  region,  solutions  of 
equations  of  the  type  (13.2)  and  (13.2)  are 

2  3 

,  -r  /l 

imfiB  -  /r“l(l-e  )gQdVa  ,  (14.2) 

‘w-=7Q:p0)av0,(i4.3) 


where 

r'-ve2+  >is+  ca . 

Now 

lim  r  -  r  , 

lxm  r*  »  0  .  (l4.4) 

V  -*•  0  1 

v  ->  0 

0 

0 

Hence, 

for  the  concentrated  force. 

(14.2)  reduces  to 

Wb  -  r-i(l-e  *)£  .  (14.5) 


In  (14.3)  the  term  in  the  integrand  of  the  form  x)%q'Eq  *8 
transformed  according  to 
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The  surface  integral  in  (l4.6)  vanishes  because  F  =  0  outside 

V  .  Also 
o 

^  /PQ.V(ri)dvQ  =  -P.„W(r).  (14.7) 

o 

Hence,  for  the  concentrated  force, 

i  *^1 

"5B0  "  l)J  .  (14.8) 

Equations  (14.5)  and  (14.8)  constitute  the  solution  of  (13-1) 

for  the  concentrated  force.  If  l2  =*  0  then  B  =  0  and 

1  c 

the  solution  reduces  to  that  found  in  [9].  If,  in  addition, 

-t2  *  0,  the  solution  reduces  to  Kelvin's  [4,  p.  183]. 

2 

15.  Surface  effects 

It  has  been  noted,  in  Section  8,  that  the  dispersion 
relations,  for  the  material  with  microstructure,  have  forms 
similar  to  those  for  a  second  order  theory  of  extensional 
vibrations  of  plates  [11].  Interesting  types  of  motion  which 
have  been  found  to  occur  at  a  free  edge  of  a  plate,  may  there¬ 
fore  be  expected  to  appear  at  a  free  surface  of  a  material 
with  micro- structure.  Examples  are:  trapping  of  energy,  at 
a  boundary,  upon  reflection  of  waves  [24];  additional  modes 
of  vibration  with  deformation  localized  at  the  boundary  [24]; 
additional  types  of  surface  waves  at  high  frequencies  [25], 


-61- 


•  Toupin  and  Gazis  [26]  have  shown  that,  according 
to  Toupin* s  generalization  of  couple-stress  theory  [8,  Sec¬ 
tion  73 >  a  free  surface,  of  a  non-centrosymmetric  material  in 
a  state  of  homogeneous  initial  stress,  will  draw  in  or  push  out 
in  a  thin  boundary  layer  as  observed  experimentally  by  Germer, 
MacRae  and  Hartman  [27].  The  equations  of  the  elastic  material 
with  micro- structure  admit  this  phenomenon  even  in  centrosym- 
metric  materials.  To  see  that  this  will  occur,  first  add  the  ini¬ 
tial  stress  terms 


c  € .  .  +  b  "V  .  . 
o  11  o'n 


to  the  potential  energy-density  (5-5)  of  the  centrosymmetric, 
isotropic  material.  If  cq  and  bQ  are  taken  to  be  constants 
of  the  same  magnitude  but  opposite  sign,  there  will  be  a  homo¬ 
geneous,  isotropic  initial  stress,  and  o^..,  with  no 

traction  t^,  according  to  (4.2)^,  across  any  surface  of  an 
infinite  body  under  no  external  forces.  However,  from  the 
equilibrium  equations  (4.1)s  (with  and  zero)  there 
must  be  non-vanishing  Hence  the  surfaces,  say  +  h, 

of  a  plate  will  have  double  tractions  acting  on  them, 

according  to  (4.2)^.  The  removal  of  these  double  tractions, 
so  as  to  free  the  surfaces,  will  result  in  deformations  governed 
by  (8.4)  with  the  right  hand  sides  set  eqqal  to  zero.  Mow, 
the  system  of  Eqs.  (8.4)  is  equivalent  to  a  sixth  order  equa¬ 
tion  of  the  form 
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l-?2  —  j 

1  dxsj\  71  dx: 


-  6. 

dx2  • 

i 


(15.1) 


where  £2  and  are  the  non-zero  roots,  probably  conjugate 

complex,  of  (8.9)  with  'j>  =  0.  Deformations  governed  by  an 
equation  of  this  form  are  localized  at  the  surface.  If  the 
plate  is  replaced  by  a  sphere,  the  analogous  solution  is  that 
for  surface  tension  on  the  sphere. 

In  the  case  of  problems  of  equilibrium  and  propagation 
of  cracks  [28],  the  equations  of  the  elastic  material  with 
microstructure  will  supply  an  energy  accompanying  the  formation 
of  new  surfaces. 

In  stress  concentration  problems,  similar  additional 
boundary  layer  phenomena  will  occur  -  over  and  above  those  by 
which  solutions  in  couple-stress  theory  [9,  10]  differ  from 
those  of  classical  elasticity.  The  added  effects  may  alter 
the  conclusions  that  were  reached  on  the  basis  of  couple-stress 
theory. 


Acknowledgement  I  wish  to  thank  Dr.  R.  A.  Toupin  for  many 
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